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ABSTRACT 

Glycolaldehyde is a simple monosaccharide sugar linked to prebiotic chemistry. Recently it 
was detected in a molecular core in the star-forming region G31.41+0.31 at a reasonably high 
abundance. We investigate the formation of glycolaldehyde at 10 K to determine whether it can 
form efficiently under typical dense core conditions. Using an astrochemical model, we test five 
different reaction mechanisms that have been proposed in the astrophysical literature, finding 
that a gas-phase formation route is unlikely. Of the grain-surface formation routes, only two are 
efficient enough at very low temperatures to produce sufficient glycolaldehyde to match the obser- 
vational estimates, with the mechanism culminating in CH3OH + HCO being favoured. However, 
when we consider the feasibility of these mechanisms from a reaction chemistry perspective, the 
second grain-surface route looks more promising, H3CO + HCO. 
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The chemistry of dense molecular cores - 
the birth sites of massive stars - is demonstra- 
bly complex, in that large molecules composed 
of several functional groups are observed to be 
present. Of particul ar interest for their astrobio- 
logica l implications ([Remiian et al.l |2004[ ISnvder 
2006) are the isomers of composition C2H4O2, 
viz. methyl formate, acetic acid, and glycolalde- 
hyde. Glycolaldehyde (CH 2 OHCHO), a simple 
monosaccharide sugar linked with the formation of 
RNA and amino acids in terrestrial en vironments 
(|Collins fc Ferrierl Il995t IWeberl Il998l ). was de- 



The mechanism of glycolaldehyde formation in 
these environments is uncertain, although it is 
becoming increasingly clear that the site of the 
formation of large organic molecu les is the icy 



surfa ces of astronomical dust (e.g., iGarrod et al 



2006). As suggested by early models of grain- 



surface chemistry, much of the development of 
complex molecules is through fairly rapid hy- 
droge nation of frozen-out ga s-phase molecules 



tected first towards the G alactic Centre mol ecular 
cloud Sagittarius B2(N) (|Hollis et al.ll2000l) . and 
more recently towards a star- forming hot m olecu- 
lar core, G31.41+0.31 (|Beltran et al.ll2009l ). both 
rich sources of molecules. 



(e.g., iTielens fc WhittetJ Il997l) . Once the end- 
points of these processes have been reached (e.g., 
C— >CH 3 OH, N-J>NH 3 ) molecules and radicals 
must move through the ice lattices in order to 
build the large organic molecules we detect in re- 
gions of star formation. 

In this paper we investigate the formation of 
glycolaldehyde in a collapsing cloud core at 10 K 
by comparing five mechanisms that have been sug- 
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Table 1 

Summary of proposed reaction pathways 





Reaction 


Reference 


Medium 


Method 


Al. 


g-H 2 + Yw > g-OH + g-H 








A2. 


g-CH 4 + — > g-CH 3 + ff-H 








A3. 


g-CH 3 + g-OH — > g-CH 3 OH 


Sorrell (2001) 


grain mantle 


theory 


A4. 


g-CO + g-H — ► g-HCO 




(H 2 0/CH 4 /NH 3 /CO) 




A5. 


g — > . 1 1 : ; v - 1 1 1 -f- g — 1 1 v. 1 ' ^ g — v. 1 1 ; 1 / 1 1 1. 1 1 1 / -f- y — 1 1 








Bl. 


g-CH 3 OH + CRP — ► g-CH 2 OH + g-H 








B2. 


g-CO + g-H — ► g-HCO 


Bennett & Kaiser ( 2007b) 


grain mantle 


experiment 


B3. 


g-CH 2 OH + g-HCO — > g-CH 2 OHCHO 




(CH3OH/CO) 




CI. 


H+ + H 2 CO — > H 2 COH+ + H 2 








C2. 


H 2 COH+ + H 2 CO — > CH 2 0HCH 2 0+ 


Halfen et al. (2006) 


gas 


theory 


C3. 


CH 2 OHCH 2 0+ — > CH 2 OHCHOH+ 








C4. 


CH 2 OHCHOH + — ► CH 2 OHCHO + H+ 








Dl. 


g-CO + g-H + g-H — ► g-H 2 CO 








D2. 


g-CO + g-H — ► g-HCO 


Beltran et al. (2009) 


surface 


theory 


D3. 


g-H 2 CO + g-HCO + g-H — ► g-CH 2 OHCHO 








El. 


g-CO + g-H — ► g-HCO 








E2. 


g-HCO + g-C — ► g-HC 2 








E3. 


g-H0 2 O + g-H — ► g-CH 2 CO 


Charnlev & Rodsers (2005) 


surface 


theory 


E4. 


g-CH 2 CO + g-H — > g-CH 2 CHO 






E5. 


g-CH 2 CHO + g-O — ► g-OCH 2 CHO 








E6. 


g-OCH 2 CO + g-H — > g-CH 2 OHCHO 









Note. — g— signifies a grain-surface species, \\v signifies a UV photon and CRP signifies a cosmic ray particle. 



gested in the astrophysical literature. These mech- 
anisms are highly speculative; all are without as- 
sociated reaction rate coefficients and many of the 
reactions involved have not previously been in- 
cluded in astrochemical models to assess their ef- 
fectiveness. We aim to constrain the possible for- 
mation routes of glycolaldehyde in cold cores, by 
investigating the wide parameter space resulting 
from the lack of existing constraints. This work is 
an initial investigation which forms part of a larger 
program looking into the formation of glycolalde- 
hyde in the dense interstellar medium, using the 
combined tools of astrochemical modelling, exper- 
imental surface chemistry and quantum chemical 
calculations. 

Dense prestellar cores have a very limited range 
of temperatures, from ~ 7-1 IK (e g..|Pagani et al 



2007[lBergin et al.l2006ULai et al.f2003HHotzel et al.l 

2002) , and there is experimental evidence (|6berg et al 
2009 : iBennett fc Kaiserl l2007bl ) that glycolalde- 
hyde forms in such low temperature environments. 
Given the uncertainties in reaction rates, we look 
at a large parameter space, and conservatively 
restrict ourselves to simple hydrogenation of the 
species which are frozen out onto grain surfaces. 
In this way we identify which of the mechanisms 



suggested in an ad hoc manner are feasible for the 
production of glycolaldehyde in molecular cores 
such as G31.41+0.31. In S}2] we give details on 
the selected mechanisms which we investigate. 
gives an overview of our model, and the procedure 
which wc follow in investigating the mechanisms. 
In 2] and ij5] we draw out some results from our 
modelling, and evaluate them bearing the chemi- 
cal energetics of the reactions in mind. In Sj6] we 
conclude with a summary of our findings. 

2. Proposed pathways to glycolaldehyde 

Since its detection in space, there has been 



(e.g. Sorrelj 2001: Charnlev & Rodgers 


2005|: 


Halfen et al.l2006HBennett & Kaiserl2007b 


: Beltran et al. 


2009. and others) . Several mechanisms have been 



proposed, including both gas-phase and surface re- 
actions, and experiments have been conducted on 
laboratory surface analogs. We summarise some 
of the work which has been carried out below, and 
in Table [T] We do not consider hig h temperature 



(~300K) formation routes, e.g., iJalbout et al 
(|2007l) . 

Some of these reactions have been tested in hot 
core models, at temperatures up to 200 K. For 
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example, Garrod et al. ( 20081 ) (and, presumably, 
Laas et al.ll201ll ) incorporate reactions B3 and E6. 
However, the chemistry in these warm tempera- 
ture regimes is somewhat different, since surface 
radicals can be sufficiently energetic to overcome 
diffusion barriers, affording them greater mobility 
on grain surfaces. We only consider these reactions 
at 10 K in order to test whether glycolaldehyde can 
form efficiently in the isothermal collapse phase of 
star formation. 

Below we summarise the work from which the 
reaction mechanisms in Table [T] comes. 

2.1. Mechanism A 



2.3. Mechanism C 



Sorrelll (|200ll ) discusses the theory of processing 



icy grain mantles in the interstellar medium with 
ultraviolet (UV) radiation, producing high con- 
centrations of free radicals (particularly OH and 
CH3). These radicals then react in the grain man- 
tles in order to produce large organic molecules 
such as amino acids and sugars, with the energy 
for thermal hopping coming from grain-grain col- 
lisions. The resulting large organics (including 
glycolaldehyde) would then be desorbed into the 
gas phase following mantle explosions, and despite 
some fraction of these large molecules being de- 
stroyed in the process, some would remain intact. 

2.2. Mechanism B 



Bennett fc Kaiserl (|2007bl ) simulated the bom- 



bardment of grain mantles with cosmic ray par- 
ticles by irradiating laboratory methanol/carbon 
monoxide ices with energetic electrons at 11 K. 
Cosmic rays can penetrate entire grains, produc- 
ing up to 100 suprathermal particles each, which 



then i onise (methan e) ice molecules (jKaiser et al 
19971 : iKaLs er 2002). The resulting high-energy 
electrons (~5keV) may then affect the man- 
tle chemistry by forming radicals, which subse- 
quently react to form large organic molecules. In 
a methanol/carbon monoxide ice these large or- 
ganics include C2H4O2 isomers. The experiment 
showed that both glycolaldehyde and methyl for- 
mate were formed, in addition to many smaller 
molecules and radicals. Acetic acid was not de- 
tected, but can be formed in meth ane /carbon 
dioxide ices ( Bennett fc Kaiserll2007al) . 



Halfcn et al.l (|2006l ) postulate that glycolalde- 
hyde may be formed in the gas phase through 
acid-catalysed reactions of form aldehyde, based 



on research on formose reactions (Butlerow 1861 



Breslowlll959h . Formaldehyde would react with its 
protonated form to create an intermediate species, 
which would then undergo reorganisation into pro- 
tonated glycolaldehyde. There is some experimen- 
tal evidence for this method, although it is unclear 
whether the re sulting C2H4O2 iso mer is in fact 
glycolaldehyde (jjalbout et al.ll2007l ). 



2.4. Mechanism D 



Beltran et al. ( 20091 ) highlighted the potential 
importance of the HCO radical in glycolaldehyde 
formation. They suggested that reactions between 
HCO and methanol (or methanol derivatives) or 
formaldehyde could occur rapidly on grain sur- 
faces in hot cores. Gas-phase routes would be too 
inefficient. The simplicity of the reaction pathway, 
which is driven by rapid hydrogenation and the re- 
action of small surface radicals, means that glyco- 
laldehyde formation could be efficient when den- 
sities are high. Only small amounts of CO would 
need to be processed on grains. 

2.5. Mechanism E 



Charnlev fc Rodgersl (|2005l ) suggested that 
complex molecules build up on grain surfaces 
through the aggregation of common atoms, since 
at low temperatures only atoms are likely to be 
mobile. At early times, atoms such as C, N or O 
may accrete significantly, whereas at late times, 
when most heavy atoms will have frozen out, hy- 
drogenation of molecules will dominate. Such a 
scheme could not only lead to the formation of 
glycolaldehyde, but also to other large molecules 
such as acetic acid and aminomethanol. Methyl 
formate, however, cannot be formed through this 
kind of pathway, but only through the combina- 
tion of relatively large surface radicals. 

3. The chemical model and scientific pro- 
cedure 

In order to test which of the suggested routes to 
glycolaldehyde formation are feasible in dense core 
environments, we have incorporated the above 
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chemical reactions (Table Q} into a model of a 
hot molecular core. The model, described be- 
low in mo r e det ail, is based on that described in 
IViti et al.1 (12004 . 

3.1. The model 

The model is a two-phase time-dependent 
model which follows the collapse of a prestellar 
core (phase I), followed by the subsequent warm- 
ing and evaporation of grain mantles (phase II). 
We only consider phase I, since in this work we 
wish to investigate the potential formation of gly- 
colaldehyde at low temperatur e s, as suggested by 
experi ment ([Oberg et al. 2009; iBennett &: Kaiserl 
2007b|). In phase I, a diffuse cloud of density 
10 2 molecules cm~ 3 undergoes free-fall collapse 
until it has reached a density of ~10 7 cm~ 3 . This 
occurs on a timescale of half a million years and at 
a temperature of 10 K. During the collapse, atoms 
and molecules collide with, and freeze on to, grain 
surfaces. We assume that hydrogenation occurs 
rapidly on these surfaces, so that, for example, 
some percentage of carbon atoms accreting will 
rapidly become frozen-out methan e, CH4. Initial 
atomic abundances are taken from ISofia fc Meveil 
(|200ll) . as in IViti et all (|2004) . We_ employ the 



reaction rate data from the UDfA0(J3 astrochem- 
ical database, augmenting it with grain-surface 
(hydrogenation) reactions and those reactions in- 
cluded in Table [TJ In the formation of glyco- 
laldehyde we only consider the most propitious 
of circumstances. We discount the destruction of 
glycolaldehyde on the grains through cosmic ray 
strikes, photodissociation or further reaction, so 
that the quantity of glycolaldehyde formed can 
be regarded as an upper limit. We also assume 
that due to the cold temperature, no species will 
desorb from the grains except H and He. This as- 
sumption is reasonable when considering thermal 
desorption, but neglects the effect of non-thermal 
desorption mechanisms on both glycolaldehyde 
and the reactants which go into its formation. 
Any proposed reaction pathways that do not pro- 
duce reasonable amounts of glycolaldehyde under 
these conditions can surely be dismissed from con- 
sideration. 



3.2. The procedure 

Firstly, we investigate the five schemes individ- 
ually, varying key model parameters, such as the 
rate coefficients, the final collapse density, the in- 
cident UV field, the cosmic ray ionisation rate and 
the hydrogenation efficiency of accreted molecules. 
Finally we look at all the mechanisms together, to 
find which is the most efficient in competition. 

Many of the rate coefficients of the reactions in 
Table [T] are completely unknown. In light of this, 
we consider a wide parameter space, covering up 
to 14 orders of magnitude in reaction rate. Our 
aim is to understand the behavior of the reactions 
in each mechanism, and what effect they have on 
the abundance of glycolaldehyde, not to determine 
reaction rates. However, through our investigation 
we may be able to better constrain possible reac- 
tion rates. Where practical, we adopt identical 
or similar gas-phase reaction rates from UDfA06 
or KiDAEI for unknown grain-surface rates as a 
conservative initial estimate (Table [3]), given that 
the grain surface is thought to act as a catalyst. 
Glycolaldehyde is not included in the standard 
UDfA06 database, so we utilise rates from anal- 
ogous reactions which produce methyl formate, or 
we make very conservative estimates. In varying 
the rates, we vary only the a-parameter in the rate 
coefficients: k = a(T/300K) /3 exp(- 7 /T) for two- 
body reactions, k — a for reactions with cosmic 
rays, and k = aexp(— jAy) for photoreactions. 

We test two final collapse densities, rtf = 10 6 
and 10 7 cm~ 3 , which has implication s also for 



the fr eeze-out percentage of molecules. IViti et al 



1 http://www.udf a. net 



(|200ll ) argued that in hot cores freeze-out is never 
total and in fact some gaseous CO is always ob- 
served, even in region s where no millimetre contin- 
uum is detected (e.g.. lMolinari et al.ll2000j ). Hence 
we impose the constraint that a maximum of 90% 
of the circumstellar material is frozen out at a final 
density of nf = 10 7 cm~ 3 , and 75% at 10 6 cm~ 3 . 

The strength of the impinging UV field within 
the core was adjusted solely for mechanism A, 
since it involves the UV processing of molecules 
on grains. We investigated the effects of scaling 
the standard interstellar UV field strength, Go, by 
up to 30 times. When investigating reaction Bl, 
we increase the cosmic ray ionisation rate in the 



2 http: //kida. obs .u-bordeauxl . f r 
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Table 2 

Hydrogenation percentages for accreting species 



Accreting species 


Products in regime 1 (fl) 


Products in regime 2 (f2) 


O 


2% O, 18% OH, 80% H 2 


1%0, 9% OH, 90% H 2 


CO 


70% CO, 20% HCO, 5% H 2 CO, 5% CH 3 OH 


60% CO, 25% HCO, 10% H 2 CO, 5% CH3OH 


c 


2% C, 3% CH, 5% CH 2 , 20% CH 3 , 70% CH 4 


1% C, 4% CH, 8% CH 2 , 12% CH 3 , 75% CH 4 


HCO<+) 


20% HCO, 40% H 2 CO, 40% CH3OH 


10% HCO, 45% H 2 CO, 45% CH3OH 


OH 


10% OH, 90% H 2 


5% OH, 95% H 2 



Table 3 

Adopted reaction rates for key reactions 



Reaction a coefficient Comment 



Al. g-H 2 + hu — ► g-OH + g-H 


5.9x10" 


10 


(7= 


=1.7) 


Based on gas-phase 


A2. g-CH 4 + hv — ► g-CH 3 + g-H 


2.2x10" 


10 


(7= 


=2.2) 


Based on gas-phase 


A3. g-CH 3 + g-OH — ► g-CH 3 OH 


1.0x10" 


9 






Based on similar gas-phase 


A4. g-CO + g-H — ► g-HCO 










Assumed fast 


A5. g-CH 3 OH + g-HCO — > g-CH 2 OHCHO - 


- g-H 3.0x10" 


17 






Methyl formate rate retarded by x50 


Bl. g-CH 3 OH + CRP — ► g-CH 2 OH + g-H 


1.3x10" 


17 






Standard £ (CRP = cosmic ray) 


B2. g-CO + g-H — > g-HCO 










Assumed fast 


B3. g-CH 2 OH + g-HCO — ► g-CH 2 OHCHO 


3.0x10" 


17 






Conservative estimate 


CI. H+ + H 2 CO — ► H 2 COH+ + H 2 


6.3x10" 


9 






Based on similar gas-phase 


C2. H 2 COH+ + H 2 CO — > CH 2 OHCH 2 0+ 


1.0x10" 


10 






Based on similar gas-phase 



C3. CH 2 OHCH 2 0+ — > CH 2 OHCHOH+ 






Assumed fast 


C4. CH 2 OHCHOH+ — ► CH 2 OHCHO + H+ 






Assumed fast 


Dl. g-CO + g-H + g-H — ► g-H 2 CO 






Assumed fast 


D2. g-CO + g-H — > g-HCO 






Assumed fast 


D3. g-H 2 CO + g-HCO + g-H — > g-CH 2 OHCHO 


3.0x10" 


-17 


Conservative estimate, based on 
methyl formate reaction 


El. g-CO + g-H — ► g-HCO 






Assumed fast 


E2. g-HCO + g-C — > g-HC 2 


1.0x10" 


-10 


Based on similar gas-phase 


E3. g-HC 2 + g-H — ► g-CH 2 CO 


1.0x10" 


-9 


Based on similar gas-phase 


E4. g-CH 2 CO + g-H — > g-CH 2 CHO 


5.0x10- 


-9 


Based on similar gas-phase 


E5. g-CH 2 CHO + g-O — > g-OCH 2 CHO 


1.0x10" 


-12 


Based on H-atom gas phase reaction, 
retarded 


E6. g-OCH 2 CO + g-H — ► g-CH 2 OHCHO 


3.0x10" 


-17 


Conservative estimate, based on 
methyl formate reaction 



Note. — Reactions A4, B2, D2 and El are identical. 
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model globally by up to 1 000 times. 

Finally, we tested two different grain-surface 
hydrogenation regimes, one where the products 
were more saturated ("f2"), and one less satu- 
rated ("fl"),. These regimes are represented in 
Table [2j Hydrogenation rates on grain surfaces are 
unknown, but the calculation of hydrogen-atom 
hopping and tunnelling rates show that they are 
rapid in comparison to other grain-surfa ce reac- 
tions (jGoumans et alJl2007trTielenslll989l) . At low 
temperatures, the recombination of physisorbed 
atomic hydrogen with chemisorbed atoms domi- 
nates. At temperatures <20K, molecular hydro- 
gen for mation efficiency o n grain surfaces is near 
unity (|Cazaux fc Tielensl l2004t IWilliams et al. 
20071 ). We assume that hydrogenation of species 
on grain surfaces is instantaneous. 

We combine these free parameters into a grid 
of model results. 

4. Results 

We have calculated approximately 450 models 
to investigate the formation of glycolaldehyde at 
10 K during the isothermal collapse phase of star 
formation. Many of these permutations arise from 
varying the a parameter in the rate coefficient of 
the reactions involved in the mechanisms by a fac- 
tor of up to 10 ±7 . 

4.1. The effect of scaling the UV field 
strength on mechanism A 

In order to explore the effects of enhanced UV 
irradiation of surface ices in mechanism A (Ta- 
ble [T|), we have used the standard reaction rates 
as shown in Table [3l and the more conservative 
hydrogenation regime (Table [2]). We varied the 
strength of the UV field by up to a factor of 30 
over the standard interstellar field, and the re- 
sults are plotted in Fig. [TJ Intuitively one would 
expect that increased grain processing of surface- 
bound H2O and CH4 would lead to a greater pro- 
duction of glycolaldehyde. However, higher UV 
fluxes mean that the gas-phase species which go 
on to form CH4 in particular on the grains are 
destroyed by reaction with abundant photodisso- 
ciation products like H^ and H + . Thus the limited 
abundance of grain-surface CH4 limits the forma- 
tion of glycolaldehyde. Water ice increases very 
marginally in abundance when the core is under 



higher levels of irradiation. 

4.2. Reaction rate analysis of mechanism 
A 

Results of the model can be found in Fig. [2] 
where we plot the fractional abundance of gly- 
colaldehyde obtained at the final density of the 
collapse, n,{, against the scale factor of the "stan- 
dard" rates, as found in Table [3J We compare 
models with differing rif and in the two differ- 
ent hydrogenation regimes, fl and f2. The frac- 
tional abundance of glycolaldehyde in the model 
is highly dependent on the rate of the final reac- 
tion of mechanism A, A5 - it scales linearly with 
the rate until 1-10 times the standard rate, after 
which the curve turns over to a maximum abun- 
dance of x(CH 2 OHCHO)~10- 5 . The effect of the 
differing hydrogenation regimes is minimal. The 
fractional abundance of glycolaldehyde is fairly in- 
sensitive to the rates of reactions A1-A3, increas- 
ing by only a factor of 10 or so whilst the rate co- 
efficients span a range of 10 14 . Reactions Al and 
A2 compete with grain surface hydrogenation of 
O and C in the provision of OH and CH3 radicals, 
respectively. As core density increases, photons 
become increasingly absorbed, meaning that reac- 
tion A3 proceeds with reactants that are products 
of grain-surface hydro genation rather tha n grain- 
surface photolysis (cf. IPeeters et al.ll2006l ). Reac- 
tion A3 itself is competing with the hydrogenation 
of adsorbed CO, which dominates at high densi- 
ties, meaning that the abundance of glycolalde- 
hyde is relatively independent of reactions A1-A3 
in general. 

Reaction A3 has been studied in the litera- 
ture, as part of investigations into methanol for- 
m ation. Experiments o n H2O-CO ice at T < 20 K 
bv lHidaka et al. ( 20041 ) do not show any evidence 
of methane production, which implies that abun- 
dances of CH 3 are low. It seems likely that 
methanol production results mainly from the suc- 
cessive hydrogenation of CO (i.e., CO — > HCO 
H 2 CO — > CHcjO — > CHcjOH ) rather than 



via reaction A3 (jHidaka et al.ll2004h . 

Reaction A4 has also been well-studied in the 



litera ture (e.g-. lHudson fc Moorell999tlWatanabe et al 
2003), and is a crucial reaction in many of the 
mechanisms studied here (it is identical to B2, 
D2 and El). HCO ice has not been detected 
in the interstellar medium, implying that its for- 



G 



Mechanism A, n, = 10 s 



Mechanism A, n, = 10 7 




Fig. 2. — The production of glycolaldehyde via mechanisms A, B, D & E, for nf=10 6 and 10 7 cm -3 (see ij4.4l 
for a description of mechanism C). The solid lines show results using hydrogenation regime fl; the dashed 
lines, regime f2. Standard rates can be found in Tabled See the electronic edition of the Journal for a color 
version of this figure. 
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mation is slower than subsequent reactions, e.g., 
H + HCO — > H2CO, which have lower acti- 
vation energies (jWatanabe fc Kouchil 120021 ) . In- 
deed, the formation of H2CO i n this way has 
been shown to be barrier less (Goumans et al 



2007). CO + H has a barrier of several thou- 



sand Kelvin in th e gas phase (see summary by 



Hidaka et al.l 120071 ). but barriers on a surface de- 



pend on_thej^omjDOj^io^ 



face (jWatanabe et alj2004lGoumans et ai1l2008l) . 
with in so me cases the reaction b ei ng completely 



barrier less (jGoumans et al.l 120081 ) . iHidaka et al 



(|2007l ) calculate a rate for A4 from an experi- 
ment involving various combinations of CO and 
H2O ice. They find that the product knnn is 
~5xl0~ 3 s _1 , which for 10~ 3 H atoms per square 
centimetre of surface (nn; typical for a large grain) 
gives a reaction timescale on the order of seconds. 
For small grains, the timescale could be on the 
order of a year (i.e., fast by astrophysical stan- 
dards). 

Reaction A5 requires the diffusion of molecules 
and radicals across a grain surface, which is 
a slow process at 10 K. However, experimen- 
tal results show that there is some evidence of 
complex surface reaction s, even at 10 K (e.g., 
IWatanabe fc Kouchill2002l ) . Figure[2]suggests that 
even at rates slower than our conservative stan- 
dard rate, significant amounts of glycolaldehyde 
form via this mechanism. 

Increasing nf by an order of magnitude from 
10 6 to 10 7 cm -3 has the effect of increasing gly- 
colaldehyde production for a given reaction rate 
by a corresponding order of magnitude, approx- 
imately. This reflects the greater collisional 
rate between molecules and grains, and thus a 
greater freeze-out rate. The peak fractional abun- 
dance of glycolaldehyde produced via mechanism 
A, x(CH 2 OHCHO)-10- 5 , is unaffected by the 
change in nf , showing that a significant proportion 
of the available carbon ends up in glycolaldehyde 
at the most extreme rates investigated, something 
which is unlikely to occur naturally. 

4.3. Reaction rate analysis of mechanism 
B 



tions where the cosmic ray ionisation rate is 
increased by seven orders of magnitude, less 
than x(CH 2 OHCHO)~10- 9 results (Fig. EJ. The 
rates adopted for reactions Bl and B3 are crit- 
ical to the amount of glycolaldehyde produced, 
with the yield scaling linearly with the adopted 
rate. Using our standard rates for Bl and B3, 
z(CH 2 OHCHO)~10- 13 -10- 14 , sign ificantly lower 



The reaction mechanism suggested bv lBennett fc 
(|2007bh . which was identified experimentally in 
the laboratory, is very inefficient at produc- 
ing glycolaldehyde at 10 K. Even under condi- 



Kaiser 



than that observed in G31.41+0.31 (jBeltran et al 
2009). Moreover, we do not include the hydro- 
genation of the hydroxymethyl (CH2OH) radical 
to methanol in our reaction scheme, which surely 
must be rapid and compete with reaction B3. 

4.4. Reaction rate analysis of mechanism 
C 

Mechanism C, the only gas-phase reaction 
mechanism we investigate, also is not particu- 
larly efficient in the production of glycolalde- 
hyde, producing a;(CH 2 OHCHO)<10- 10 at the 
most enhanced values of the reaction rate coeffi- 
cients. We do not increase the rate coefficient of 
reaction CI beyond 6.3x 10~ 7 cm 3 s _1 , since this 
would be incredibly fast for a gas-phase reaction. 
In fact, we only vary the rate of CI for com- 
pleteness, since the standard reaction rate (from 
Tanner et al.l Il979h is accurate to 25% accord- 
ing to the UDfA database. Both reactions CI 
and C2 produce linearly increasing amounts of 
CH2OHCHO with increasing reaction rate coef- 
ficient, producing a;(CH 2 OHCHO)=10- 13 at the 
standard rates. The limiting factor in this mecha- 
nism is the availability of gas-phase formaldehyde, 
which at 10 K is only 1% of the total formaldehyde, 
the rest being frozen onto grain surfaces. However, 
we do not include non-thermal desorption mech- 
anisms in our simple model, which could increase 
the amount of form aldehyde in the gas phase. 
Roberts et al.l (|2007l ) show that non-thermal des- 
orption can return a large proportion of H2CO 
to the gas phase in extreme cases. It is not clear 
how their results would apply to our situation, 
where the density is larger, and thus freezeout 
more rapid. 

4.5. Reaction rate analysis of mechanism 
D 



This grain-surface re action based on the work 
of iBeltran et al.l (|2009l ) involves the products of 
rapid hydrogenation of CO and HCO^', meaning 



S 



that the only instructive reaction rate to investi- 
gate is that of reaction D3. We simplify the three- 
body reaction proposed by iBeltran et al. ( 20091 ) 
with a two-body reaction by assuming that the 
H-atom addition is rapid. A maximum fractional 
abundance of glycolaldehyde of ^10~ 6 is produced 
when the standard rate for the reaction is in- 
creased by 1000-10000, depending on nf. At the 
standard rate, x(CH 2 OHCHO)«10- 9 "- 10 , some- 
what smaller the derived from their observations, 
but within their error constraints (see Sect. [5]). 

4.6. Reaction rate analysis of mechanism 
E 

The atom-addition mec hanism suggested by 
Charnlev fc Rodgersl (|2005f) is the most complex 
that we have considered, involving six two-body 
reactions. Three of these reactions are reactions 
with H, which are assumed rapid, but the surface 
migration of heavier atoms is significantly slower 
due to the large diffusion barri ers involved (e.g., 
Leitch-Devlin fc Williamsl 119841 ) . The rates of re- 
actions E2-E4 have little effect on the abundance 
of glycolaldehyde, showing largely flat profiles in 
Fig. [2] The gradient of reaction E5, where an oxy- 
gen atom is added to the molecule, is somewhat 
steeper, but reaction E6 is the crucial reaction in 
this mechanism. We have assumed a fairly conser- 
vative value of 3xl0 _17 s _1 for this reaction, due 
to a potential barrier in the H-addition process, 
but if the reaction were 1Q 3 -- 5 times faster than 
expected, it could produce a;(CH 2 OHCHO)-10 _6 . 

4.7. Further experimentation 

The analysis performed thus far has been some- 
what artificial, since if the reactions contained in 
the suggested mechanisms occur, then they likely 
occur in competition with each other (and other 
reactions), i.e., reactions in other mechanisms are 
not 'switched off'. To investigate this, we per- 
formed two further experiments where we set the 
rates of all the reactions considered to their stan- 
dard rates, and where we set them to "optimal" 
rates. By optimal here we mean adopting a values 
for where the rate profiles in Fig. [2] turn over to 
become flat. Rates for mechanism C were kept as 
standard, since the rate profiles do not turn over. 

We find that when using standard rates for 
all reactions, x(CH 2 OHCHO)= 1.5 x 10~ 7 . In 
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Fig. 1. — The dependance of glycolaldehyde abun- 
dance upon UV radiation field intensity, for mech- 
anism A. See the electronic edition of the Journal 
for a color version of this figure. 
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Fig. 3. — The formation route of glycolaldehyde 
with all mechanisms in operation, using standard 
reaction rates. See the electronic edition of the 
Journal for a color version of this figure. 
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this competitive environment, reactions A5 and 
E6 are the dominant routes to the formation of 
glycolaldehyde, with A5 being the most efficient 
by a considerable margin as the core collapses 
(see Fig. [31 left). Using "optimal" rates an ex- 
tremely large amount of glycolaldehyde can form, 
x(CH 2 OHCHO)= 1.1 x 1(T 5 , which is approxi- 
mately 6% of the elemental carbon abundance of 
the core. In this case the most dominant reac- 
tion for its formation at later times is B3 (the rate 
of which has been enhanced 10 5 times); E6 dom- 
inates at early times, with a rate enhancement of 
10 7 . Given these extreme rate enhancements, this 
scenario is unlikely. 

5. Discussion 

Having investigated the formation of glyco- 
laldehyde via five different reaction mechanisms, 
it is clear that considerable quantities of glyco- 
laldehyde can be produced. Also, very small 
amounts of glycolaldehyde can result from mech- 
anisms which are inefficient at 10 K (e.g., mecha- 
nism B). Many of the rates involved in these mech- 
anisms are completely unknown, and others have 
a large degree of uncertainty. Given that frac- 
tional abundances of glycolaldehyde can vary over 
ten orders of magnitude or more in the models, 
further work needs to be done in order to fully 
understand how glycolaldehyde forms at low tem- 
peratures in dense molecular cores, like those in 
G31.41+0.31. We are currently undertaking Den- 
sity Function Theory calculations for these five 
mechanisms, with the results to be forthcoming 
in a future publication. 

The fact that glycolaldehyde has been detected 
in G31.41+0.31 pr ovides us with some constraints 
on our modelling. iBeltran et al. I (l2009h estimated 
that the fractional abundance of glycolaldehyde 
towards this region was on the order of 10~ 8 . Due 
to the uncertainties of temperature and column 
density in their measurements, the errors in this 
estimate are large: potentially two orders of mag- 
nitude (M. Beltran, priv. comm.). This lower 
limit of x(CH 2 OHCHO)>10- 10 effectively means 
that we can reasonably exclude mechanisms B 
and C from further consideration, since they do 
not produce enough glycolaldehyde even under the 
most favourable conditions. Mechanism E only 
produces required abundances if the rate of the 



final reaction in the scheme is enhanced. Mecha- 
nisms A and D have greater fecundity, and Fig. [3] 
shows that mechanism A is considerably more ef- 
ficient when in competition. Thus it appears that 
for the rates we have adopted, the main formation 
mechanism for glycolaldehyde at 10 K is: 



CH 3 + OH 
CO + H 
CH 3 OH + HCO 



CH3OH 
HCO 

CH 2 OHCHO 



H 



on grain surfaces. The formation of both CH3OH 
and HCO in ices has been we l l-stud i ed experimen- 
tally 



(e.g.. IWatanabe et al l l20ol Irlidaka et al 
2004t IWatanabe et al.ll2003t IWatanabe fc Kouchi 



2002), and is related, with HCO being a cru- 
cial part in the formation of CH3OH, which 
is the terminal molecule in the hydrogenation 
process. Only energetic processes (e.g., UV or 
cosmic-ray irradiation) will decompose CH3OH, 
since H-abstraction to form H 2 CO is negligible 
(Iffidaka et alj|2004l ). 



5.1. A chemical evaluation of mechanisms 
A, D and E 

Building on the astrophysical models, we now 
consider reactions A, D and E from a physico- 
chemical perspective, taking account of the intrin- 
sic thermodynamic stability of reagents and prod- 
ucts. A full treatment of this is underway, but 
due to the computationally-expensive nature of 
such an investigation, here we only consider gen- 
eral principles. Firstly, in mechanism A, reaction 
A5 is found to be particularly efficient at produc- 
ing glycolaldehyde: 



CH3OH + HCO 



CH2OHCHO +H. 



However, A5 is the reaction of a stable molecule 
(methanol) with a reactive radical (formyl) to give 
stable glycolaldehyde and an H monoatom. The 
H monoatom is extremely reactive and even at 
low temperatures is very mobile, hence the rate 
of the backward reaction (addition of H to gly- 
colaldehyde) could be expected to be competitive 
with that of the forward reaction. Preliminary 
ab-initio calculations conducted at the coupled- 
cluster level with single, double and triple exci- 
tations (CCSD(T)) and a triple-zeta quality basis 
set shows that in fact reaction A5 is endothermic, 
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and consequently the rate of the reaction yielding 
glycolaldehyde would be very unfavourable. 
Reaction D3, 



H 2 CO + HCO + H 



CH2OHCHO, 



is a three-body reaction which has a vanishingly 
small probability of occurring, further hindered 
by the low temperatures considered here. How- 
ever, the products could be obtained by two se- 
quential two-step reactions: first, reaction of H 
with HCO yields H2CO, a barrierless process in 
the g as phase, according to published theoretical 
work (jGoumans et al.1120071 ) , which could combine 
with another H2CO molecule to yield the glyco- 
laldehyde product. However, H 2 CO is rather sta- 
ble and hence the reaction rate for the condensa- 
tion of two H 2 CO molecules could be expected to 
be rather slow and hence improbable. A second 
possibility is the reactio n of H with H2CO which 



accor ding to past work (jWoonl l2002t iSaebo et al 



Il983l ) yields H 3 CO as the kinetic product. H 3 CO 
could react with HCO to give glycolaldehyde, and 
this reaction ought to have a low barrier since both 
H3CO and HCO are reactive radical species. 

In mechanism E, the final promising path iden- 
tified from the astrophysical models, glycolade- 
hyde is assembled from a building block of CO via 
six stepwise monatomic addition reactions. Al- 
though the constituent reactions of mechanism E 
arc chemically viable, consideration of the physi- 
cal conditions and the reaction probabilities sug- 
gest that this pathway may be an unlikely source 
of glycolaldehyde. Under the conditions of a tem- 
perature of 10 K, only monatomic H is mobile; C 
and O are static and only become mobile during 
warm-up, which implies that reactions E2 and E5 
are highly unlikely to occur. Hence these reactions 
are likely to be rate-limiting. Construction of gly- 
colaldehyde via monoatomic addition also depends 
on a well-defined consecutive set of reactions. In 
this scheme, HCO reacts with monatomic C (re- 
action E2), yet monatomic H is present at higher 
abundance and is known to react without a barrier 



(jGoumans et al.1 120071 ) with HCO, the product of 
El, to give H 2 CO rather the product of E2, HC 2 0. 
Since mechanism E is composed of six reactions 
which are expected to be limited by at least two 
of those reactions (E2 and E5), we consider this 
pathway not to be very probable or efficient. 
Of the reactions identified by the astrophysi- 



cal models, it is suggested that mechanism D is 
most probable according to chemical considera- 
tions. Detailed ab-initio calculations are under- 
way to assess the influence of substrates on the 
reaction barriers and hence quantify the efficiency 
of mechanisms A, D and E. It should be noted 
that mechanisms B and C are viable from a chem- 
ical standpoint and these will also be considered 
in comparison to A, D and E to assess the most 
viable scheme. 

6. Summary 

We have investigated five reaction pathways for 
the formation of glycolaldehyde, a simple sugar, 
which have been previously suggested in the as- 
trophysical literature, but not thoroughly tested 
or justified. By means of a chemical model of 
an isothermally-collapsing molecular core, we have 
determined that under the physical conditions 
assumed, mechanisms B and C (suggested by 
iBennett fc Kaiserl l2007bt lHalfen et alTl2006l re- 
spectively) are relatively inefficient and are un- 
likely to be the major pathways to the formation of 
gl ycolaldehyde. Mechanisms D and E (suggested 
bv lBeltran et al.ll2009t Icharnlev fc Rodgersll2005l 
respectively) can reach observed fractional abun- 
dances of glycolaldehyde if reaction rates are en- 
hanced by factors of 100 or more over those we 
have chosen as s tandard . Fina lly, mechanism A, 
from the work of Sorrelll ( 2001 ). can produce gly- 
colaldehyde very efficiently; however, in the high 
density regions around forming protostars, pho- 
tons are unlikely to penetrate enough to initiate 
the reaction scheme as it is. Instead, the initial re- 
actants, OH and CH3, are amply supplied by the 
freeze-out and subsequent hydrogenation of atoms 
and smaller radicals. 

Further evaluation of the reaction schemes tak- 
ing reagent and product stability into account 
leads us to expect that mechanism A may not be 
as efficient as the astrophysical modelling suggests 
because of a high barrier for the final step in the 
formation of glycoladehyde. However, mechanism 
D may be more likely, provided the final reaction 
in the scheme proceeds as two two-body reactions 
involving H3CO and HCO radicals, rather than 
a single three-body reaction. Detailed examina- 
tion of the reaction barriers according to high level 
quantum chemical methods is currently being un- 
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dertaken, taking into account the role of different 
substrates. 

Finally, the list of reactions considered here is 
not exhaustive and we are seeking to use astro- 
physical modelling, quantum chemical and exper- 
imental approaches to identify whether other, as 
yet unreported, pathways are more efficient at pro- 
ducing glycolaldehyde. 

The authors appreciate fruitful discussions with 
T. P. M. Goumans, which have benefitted the con- 
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